Electronic Raman scattering in non-centrosymmetric superconductors 
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We formulate a theory for the polarization-dependence of the electronic (pair-breaking) Raman 
response for the recently discovered non-centrosymmetric superconductors in the clean limit at zero 
temperature. Possible applications include the systems CePtaSi and Li2Pd I Pt3_ a; B which reflect 
the two important classes of the involved spin-orbit coupling. We provide analytical expressions for 
the Raman vertices for these two classes and calculate the polarization dependence of the electronic 
spectra. We predict a two-peak structure and different power laws with respect to the unknown 
relative magnitude of the singlet and triplet contributions to the superconducting order parameter, 
revealing a large variety of characteristic fingerprints of the underlying condensate. 

PACS numbers: 74.70.-b 74.25. Gz 74.20.Rp 



The order parameter of conventional and unconven- 
tional superconductors is usually classified as either spin 
singlet (even parity) or spin triplet (odd parity) by the 
Pauli exclusion principle. A necessary prerequisite for 
such a classification is, however, the existence of an inver- 
sion center. Something of a stir has been caused by the 
discovery of the bulk superconductor CePt 3 Si without 
inversion symmetry [lj, which initiated extensive theo- 
retical 0, Q and experimental [j, [B| studies. In such 
systems the existence of an antisymmetric potential gra- 
dient causes a parity-breaking antisymmetric spin-orbit 
coupling (ASOC) that leads to a splitting of the Fermi 
surface and, moreover, gives rise to the unique possibil- 
ity of having admixtures of spin-singlet and spin-triplet 
pairing states. At present, however, the relative magni- 
tude of both contributions to the superconducting order 
parameter is unknown. 

In this letter, we propose that inelastic (Raman) light 
scattering provides a powerful tool to solve this problem 
and, in general, to investigate the underlying condensate 
in such parity-violated, non-centrosymmetric supercon- 
ductors (NCS). This is because various choices of the pho- 
ton polarization with respect to the location of the nodes 
on the Fermi surface allow one to draw conclusions about 
the node topology and hence the pairing symmetry. An 
example for the success of such an analysis is the work 
by Devereaux et al. [Q\ in which the ci c 2_ y 2-synimetry of 
the order parameter in cuprate superconductors could be 
directly traced back to the frequency-dependence of the 
electronic Raman spectra, measuring the pair-breaking 
effect. Therefore, our predictions of the polarization de- 
pendence of Raman spectra enable one to draw conclu- 
sions about the internal structure of the parity-mixed 
condensate in a given NCS. 

The model Hamiltonian for nonintcracting electrons in 
a non-centrosymmetric crystal reads Q 





FIG. 1: (color online) The angular dependence of |gkj for the 
point groups d v and 0(423). Since dk]|gk, these plots show 
also the magnitude of the gap function in the pure triplet case 
for both point groups. 



where £k represents the bare band dispersion, a, a' =t, j, 
label the spin state and t are the Pauli matrices. The sec- 
ond term describes an ASOC with a coupling gk [a]- In 
NCS two important classes of ASOCs are realized reflect- 
ing the underlying point group Q of the crystal. This is 
illustrated in Fig. [1] Therefore, we shall be interested in 
the tetragonal point group C^ v (application to CePtaSi, 
for example) and the cubic point group (3(432) (appli- 
cable to the system Li 2 Pd x Pt3_ a; B). For Q = C± v the 
ASOC reads 01 



gk = g±(k x e z ) + g\\k x k y k z (kl - fc 2 )e 2 
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In the purely two-dimensional case (g|| = 0) one recovers 
what is known as the Rashba interaction M, fl~oi ] . For the 
cubic point group Q = 0(432) gk reads [11| 

gk = gik (3) 

- ga [k(k 2 y + k 2 z )e x + k y (k 2 z + k 2 x )e y + k z (k 2 x + k 2 y )e 

where the ratio g3/gi — 3/2 is estimated by Ref. fill ]. 

Diagonalizing the Hamiltonian one finds the eigenval- 
ues £k± = £k ± |gk|, which physically correspond to the 
lifting of the Kramers degeneracy between the two spin 
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FIG. 2: (color online) Raman spectra for a pure triplet order 
parameter (if) = 0) for Bi,2 polarization of the point group 
Civ in backscattering geometry. The ABM (axial) state with 
| d.k ] = do sin 8 is displayed in blue and the polar state with 
|dk] = do \ cos 6 1 in green. For a comparison, also the Raman 
response for the BW state (red) with |dk| = do is shown. 



states at a given k. Sigrist and co-workers have shown 
that the presence of the ASOC generally allows for an 
admixture between a spin-triplet order parameter and 
a spin-singlet pairing gap This implies that we can 
write down the following ansatz for the energy gap ma- 
trix in spin space A k(T(T ' = [(i/j k (T)l + d k (T) • r)ir v ] <T<T > > 
where tpk(T) and dk(T) reflect the singlet and triplet 
part of the pair potential, respectively. It is then easy to 
see that the ASOC is not destructive for triplet pairing 
if one assumes dk||gk [ILEl- This results in the following 
ansatz for the gap function on both bands (+,— ) [121 ]: 



A k± = tp ± d\g k \ =i>(l± p|g k |) = A± 



(4) 



where the parameter p = d/tp represents the unknown 
triplet-singlet ratio. Note that for I^Pd-jPta^B this 
parameter seems to be directly related to the substitu- 
tion of platinum by palladium, since the larger spin-orbit 
coupling of the heavier heavier platinum is expected to 
enhance the triplet contribution [13j |. This seems to be 
confirmed by penetration depth experiments [ill ]. 

The T=0 electronic Raman response in a single band 
is given by the imaginary part of (1 <-» 7k = 1) [141 ] 



X77M 



4?M 



XirV) 



(5) 



where the index 7 = 7k denotes the momentum- 
dependent Raman vertex that describes the coupling of 
polarized light to the sample. Note that the second term 
in Eq. ([5]) is often referred to as the screening con- 
tribution that originates from gauge invariance. Since 
the ASOC leads to a splitting of the Fermi surface, the 
total Raman response is given by X7° tal = Sa= 



with x^j 



over the spin variable a is replaced by a summation over 
the pseudo-spin (band) index A. With Eq. ([4]), the un- 
screened Raman response for both bands can be written 

as 



cv Y (o)± = gjVjfy 

°*77 
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where we allow for a different density of states Np 
on both bands and (. . .)fs denotes an average over 
the Fermi surface. For small momentum transfers and 
nonresonant scattering, the Raman tensor is given by 
7k = ™E a ,^ e a (d 2 e(k)/dk a dk/3)e I fj where e 5,/ denote 
the scattered and incident polarization light vectors. The 
light polarization selects elements of this Raman tensor, 
where 7 k can be decomposed into its symmetry compo- 
nents and, after a lengthy calculation, expanded into a 
set of basis functions on a spherical Fermi surface. Our 
results for the tetragonal group Q4.V 3X6 
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X 77 (A±), in which the usual summation 



-f k ,i cos sin 2fc 6 

k=a 1=0 

E 7M cos(4/ - 2)0 sin 2fc 

k=l 1=1 

)< k + l 

J2 7k,l sin(4/ - 2)0 sin 2fe 9 

k=l 1=1 

and for the cubic group 0(432) we obtain 

00 'i 2 
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7e(d = 7o(2 - 3sin 2 8) + ... 

7 S( 2) = ^ ^ 7fc,z cos(4Z - 2)0 sin 2fc 9 
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in a backscattering-geometry experiment (zz) [l5j . In 
what follows, we neglect higher harmonics and thus use 
only the leading term in the expansions of 7 k . 

Before studying the mixed-parity case, it is instruc- 
tive to analyze the Raman response for pure triplet p- 
wave pairing, see Fig. [2] Some representative examples 
are the Balian-Werthamer (BW) state, the Anderson- 
Brinkman-Morel (ABM or axial) state, and the polar 
state. The simple pseudoisotropic BW state with dk = 
dok [equivalent to Eq. |3]) for g 3 = 0], as well as pre- 
vious work on triplet superconductors, restricted on a 
(cylindrical) 2D Fermi surface, generates the same Ra- 
man response as an s-wave superconductor jl6| . In 3D 
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we obtain more interesting results for the axial state with 
dk = do(k y e x ~ k x e y ) [equivalent to Eq. ([2]) for gn = 0]. 
The Raman response for this axial state in Bx and B 2 
polarizations for Q = C^ v is given then by 
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with the reduced frequency x — lu/2cIq. An expan- 
sion for low frequencies reveals a characteristic expo- 
nent [xbj 2 oc (uj/2dof} which is due to the overlap be- 
tween the gap and the vertex function. Moreover, we 
calculate the Raman response for the polar state with 
dk = dak z e x , where one equatorial line node crosses the 
Fermi surface: 



with the trivial low frequency expansion x'b 2 ^ w/2c2q. 
Whereas the pair-breaking peaks for the BW and ABM 
state were both located at ui = 2do (similar to the Bi g 
polarization in the singlet <i-wavc case, which is peaked 
at 2A ), for the polar state this peak is significantly 
shifted to lower frequencies. The maximum is given by 
uj = 1.38g?o, which looks similar to the response for B 2g 
polarization in singlet d-wave superconductors, where it 
is also peaked atw< 2A . 

In general, due to the mixing of a singlet and a triplet 
component to the superconducting condensate, one ex- 
pects a two-peak structure in parity-violated NCS, re- 
flecting both pair-breaking peaks for the linear combi- 
nation [see Eq. Q] of the singlet order parameter ifa 
(extensively discussed in Ref. [17|) and the triplet order 
parameter dk (shown in Fig. [2]), respectively. The ratio 
p = d/ip, however, is unknown for both types of ASOCs. 

In Fig. [3]we show the calculated Raman response using 
Eq. © with g|| = 0. This Rashba-type of ASOC splits 
the Fermi surface into two bands; while on the one band 
the gap function is Ak = ip (1 +p|gk|) = A + , it is A_ = 
-0(1— p|gk|) on the other band. Thus, depending on the 
ratio p = d/ip, four different cases (see polar diagrams in 
the insets) have to be considered: (a) no nodes; (b) one 
(equatorial) line node (A„ band); (c) two line nodes (A_ 
band); and (d) two point nodes on both bands. Since the 
Raman intensity in NCS is proportional to the imaginary 
part of x*^ al = x 77 (A_) + x 77 (A + ), it is interesting to 
display both contributions separately (blue and red, re- 
spectively). Even though (except for ip = 0) we always 
find two pair-breaking peaks at oj /2ip — |1 ±p\ we stress 
that our results for NCS are not just a superposition of 
a singlet and a triplet spectra. This is clearly demon- 
strated in Fig[3]Ja), for example, in which we show the 
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FIG. 3: (color online) Theoretical prediction of the Raman 
spectra x" 7 (A_) [blue] and x" 7 (A+) [red] for Ai (solid lines) 
and for Bj.,2 (dashed lines) polarizations for the point group 
Ci v . We obtain the same spectra for the Bi and B2 symme- 
try. The polar diagrams in the insets demonstrate the four 
qualitative different cases. 



results for a small triplet contribution (p = 1/2). For 
X 77 (A_) we find a threshold behavior with an adjacent 
maximum value of x^ i2 (A_) = 7 2 7r 2 /8y / p^— 1 
and for x 77 (A + ) a zero Raman signal to twice the singlet 
contribution followed by a smooth increase and a singu- 
larity 18]. For the special case (b), where the singlet con- 
tribution equals the triplet one (p = 1), the gap function 
A_ displays an equatorial line node without sign change. 
Because of this nodal structure and strong weight from 
the vertex function (oc sin 2 9), many low energy quasipar- 
ticles can be excited, which leads to this square-root-like 
increase in the Raman intensity. In this special case the 
pair-breaking peak is located very close to elastic scat- 
tering (ui — 0.24-0). In Fig. [3](c) the gap function A_ 
displays two circular line nodes. The corresponding Ra- 
man response for p > 1 shows two singularities with dif- 
ferent low frequency power laws [x'b^ 2 (^~) ^ Cc, /^V' an d 
x" Bl 2 (A+) oc {ui/2^)— I)" 5 "]. Finally, forp ^> 1 one recov- 
ers the pure triplet cases (d) which is given analytically 
by Eq. ©. 

The Raman response for the point group 0(432), us- 
ing Eq. is shown in Fig. 2J We again consider 
four different cases: (a) no nodes; (b) six point nodes 
(A_ band); (c) six connected line nodes (A_ band); and 
(d) 8 point nodes (both bands) as illustrated in the in- 
sets. Obviously, the pronounced angular dependence of 
|gk| leads to a strong polarization dependence. Thus we 
get different peak positions for the E and T 2 polariza- 
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FIG. 4: (color online) Theoretical prediction of the Raman 
spectra x-ry(A-) [blue] and x~ry{A+) [red] for E (solid lines), 
T2 (dashed lines) and Ai [dotted line, only in (d)] polariza- 
tions for the point group 0(432). The insets display the point 
and line nodes of the gap function A_. 

tions in x" 7 (A + ). As a further consequence, the Raman 
spectra reveals up to two kinks on each band (+,— ) at 
uj/2ip = |1 ±p/4| and = |1 ±p\ 0. Furthermore, 

no singularities are present. Nevertheless, the main fea- 
ture, namely the two-peak structure, is still present and 
one can directly deduce the value of p from the peak and 
kink positions. Finally, for p 3> 1 one recovers the pure 
triplet case (d) , in which the unscreened Raman response 
is given by 

„ , v 2d / 2 \ 

X-^iu) oc — U{ 7 k . = ) 

Clearly, only the area on the Fermi surface with 
u>/2d> |gk| contributes to the Raman intensity. Since 
|gk| S [0,1] has a saddle point at |gk| = 1/4, we 
find kinks at characteristic frequencies w/2d = 1/4 and 
u>/ 2d = 1. In contrast to the Rashba-type ASOC, we 
find a characteristic low energy expansion oc (u>/2d) 2 for 
both the Ai and E symmetry, while oc (uj/2d) A for the 
T 2 symmetry [H. 

In summary, we have calculated for the first time 
the electronic (pair-breaking) Raman response in the 
newly discovered NCS such as CePtaSi (Q = Ca v ) and 
Li 2 Pd x Pt 3 _ x B (Q = 0(432)). Taking the pronounced 
ASOC into account, we provide various analytical results 
for the Raman response function and cover all relevant 



cases from weak to strong triplet-singlet ratio p. Our 
theoretical predictions can be used to analyze the under- 
lying condensate in parity-violated NCS and allow the 
determination of p. 

We thank P. M. R. Brydon and M. Sigrist for helpful 
discussions. 
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